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ABSTRACT 

The extratropical eddy momentum flux (EMF) is controlled by generation, propagation, and dissipation of 
large-scale eddies and is concentrated in Earth’s upper troposphere. An idealized GCM is used to investigate 
how this EMF structure arises. In simulations in which the poles are heated more strongly than the equator, 
EMF is concentrated near the surface, demonstrating that surface drag generally is not responsible for the 
upper-tropospheric EMF concentration. Although Earth’s upper troposphere favors linear wave propagation, 
quasi-linear simulations in which nonlinear eddy-eddy interactions are suppressed demonstrate that this is 
likewise not primarily responsible for the upper-tropospheric EMF concentration. The quasi-linear simula¬ 
tions reveal the essential role of nonlinear eddy-eddy interactions in the surf zone in the upper troposphere, 
where wave activity absorption away from the baroclinic generation regions occurs through the nonlinear 
generation of small scales. In Earth-like atmospheres, wave activity that is generated in the lower troposphere 
propagates upward, then turns meridionally, eventually to be absorbed nonlinearly in the upper troposphere. 
The level at which the wave activity begins to propagate meridionally appears to be set by the typical height 
reached by baroclinic eddies. This can coincide with the tropopause height but also can lie below it if con¬ 
vection controls the tropopause height. In the latter case, EMF is maximal well below the tropopause. The 
simulations suggest that EMF is concentrated in Earth’s upper troposphere because typical baroclinic eddies 
reach the tropopause. 


1. Introduction 

Large-scale baroclinic eddies shape the general circu¬ 
lation of Earth’s atmosphere. They are generated in mid¬ 
latitudes through baroclinic instability, propagate merid¬ 
ionally, and dissipate near their critical lines on the flanks 
of the jet streams (Randel and Held 1991). Meridionally 
propagating eddies transport (angular) momentum toward 
their generation region (Held 1975, 2000). Consequently, 
generation of large-scale eddies in midlatitudes and dissi¬ 
pation at lower and higher latitudes leads to a meridional 
momentum flux, with convergence in midlatitudes and di¬ 
vergence in the subtropics and, to a lesser extent, in polar 
regions (Fig. la). 

The eddy momentum flux (EMF) controls the structure 
of the mean zonal surface wind and of meridional cells. To 
first order in Rossby number, surface friction balances the 
EMF divergence averaged over an atmospheric column in 
the extratropics. This balance controls the strength and di¬ 
rection of the mean zonal surface winds. In the upper tro¬ 
posphere, EMF divergence is locally balanced by the Cori- 
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olis torque acting on the meridional wind, which accounts 
for the mass flux in the Ferrel cell, in the high-latitude po¬ 
lar cell, and to some extent in the tropical Hadley cell (see 
Schneider 2006 for a review). Hence, the structure of the 
EMF is fundamental to the mean state of Earth’s atmo¬ 
sphere. 

It is well known that EMF is concentrated in the upper 
troposphere, just below the tropopause (Fig. la). Eddy ki¬ 
netic energy (EKE oc + v'^) is also maximal in the upper 
troposphere (Fig. Ib). However, explaining the structures 
of EMF and EKE is not equivalent because the correlation 
coefficient between u' and v' varies spatially. The absolute 
value of the correlation coefficient between u' and v' gen¬ 
erally increases with height within the troposphere, from 
values smaller than 0.2 above the planetary boundary layer 
to values around 0.5 in wave-breaking regions near the 
tropopause. Hence, the concentration of EMF in the upper 
troposphere is stronger than that of EKE. The structures 
of EKE and EMF are not the same because EMF arises 
from the irreversible processes of eddy generation and dis¬ 
sipation, whereas EKE is nonzero even in reversible wave 
dynamics. 

The upper-tropospheric concentration of EMF is one of 
the most conspicuous features of atmospheric eddy fields. 
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Yet proposals of how it arises are scant. It has been pro¬ 
posed that friction acting on the eddies plays a role in re¬ 
ducing eddy amplitudes and meridional propagation near 
the surface, thus leading to reduced EMF near the surface 
(Held 2000; Vallis 2006). But observations of Jupiter and 
Saturn show that friction is unlikely to be important. On 
Jupiter and Saturn, EMF has been observed by tracking 
visible clouds in the upper troposphere, with EMF conver¬ 
gence in prograde (westerly) jets and divergence in retro¬ 
grade (easterly) jets (Salyk et al. 2006; Del Genio et al. 
2007). Although we do not have direct observations be¬ 
low the visible clouds, one can infer that EMF must be 
concentrated in a relatively shallow layer in the upper tro¬ 
posphere, because otherwise the implied transfer of EKE 
to mean-flow kinetic energy would, implausibly, exceed 
the total energy available to drive the flow (Schneider and 
Liu 2009; Liu and Schneider 2010). Indeed, simulations of 
Jupiter’s and Saturn’s atmospheres exhibit EMF concen¬ 
trated in the upper troposphere, near the tropopause, much 
like in Earth’s atmosphere, although in the simulations this 
is far above any drag layer, as it likely is for the actual giant 
planets where drag may arise magnetohydrodynamically 
in the planet’s interior (Schneider and Liu 2009; Liu and 
Schneider 2010). So upper-tropospheric concentration of 
EMF is not unique to Earth’s atmosphere but appears to 
be an ubiquitous feature of planetary atmospheres. And 
friction is unlikely to be generally responsible for the con¬ 
centration, as it does not seem to play a role in giant planet 
atmospheres. 

Understanding the EMF concentration requires under¬ 
standing the generation, propagation, and dissipation of 
wave activity, as pioneered in studies of baroclinic life- 
cyles, which revealed the central role of baroclinic growth 
followed by barotropic decay (Simmons and Hoskins 
1978, 1980; Held and Hoskins 1985; Thorncroft et al. 
1991). Wave activity generation, propagation and dissipa¬ 
tion can be diagnosed using cross sections of the Eliassen- 
Palm (EP) flux (Edmon Jr. et al. 1980). For dynamics 
that are locally quasigeostrophic (QG), the meridional and 
vertical components of the zonal-mean EP flux in pressure 
coordinates take the form 

/ \ 

F = Rcosij) I I . (1) 

V fVW/dpB J 

Here, R is Earth’s radius, (j) is latitude, / the Coriolis pa¬ 
rameter, 6 potential temperature, and p pressure. Primes 
denote departures from the zonal average (•). The wave 
activity A obeys 

-^+W-F = ^, (2) 

where ^ includes all non-conservative terms. Equation 
(2) has been shown to hold for small-amplitude waves. 


in which case the wave activity A equals the pseudomo¬ 
mentum (Andrews and McIntyre 1976), and for flnite- 
amplitude waves (Nakamura and Zhu 2010). When WKBJ 
theory is applicable, F c^A is the advective flux of wave 
activity that is carried by the group velocity vector of 
the waves (Lighthill and Lighthill 1960; Hayes 1977). In 
general, F indicates wave activity propagation. In a sta¬ 
tistically stationary state, divergence of F indicates wave 
activity generation, and convergence of F wave activity 
dissipation. 

The meridional component of the EP flux is the merid¬ 
ional eddy angular momentum flux, which is the topic of 
this paper. In this sense, EMF is linked to meridional 
wave activity propagation. EMF divergence implies con¬ 
vergence of the meridional wave activity flux, or wave 
activity dissipation. The EP flux provides a framework 
for connecting EMF divergence in the upper troposphere 
to the upward propagation of wave activity generated at 
lower levels. 

It has been suggested that EMF is concentrated in the 
upper troposphere because potential vorticity (PV) gradi¬ 
ents there are greater than in the lower troposphere, which 
favors meridional wave propagation (Held 2000, 2007). 
The larger PV gradients and stronger zonal winds aloft 
lead to a wider region in which Rossby wave refractive 
indices are positive and meridional wave propagation is 
possible. That is, the critical latitudes, near which merid¬ 
ionally propagating wave activity dissipates (Randel and 
Held 1991), are farther away from the wave activity gen¬ 
eration region (Thorncroft et al. 1991). The weaker PV 
gradients and weaker zonal winds close to the surface, it is 
suggested, lead to nonlinear saturation of baroclinic eddies 
close to their generation latitude, precluding substantial 
EMF across latitudes. But the increase with height of the 
width of the region that allows meridional wave propaga¬ 
tion is more gradual than the peaked structure of EMF: the 
meridional distance between the critical latitudes widens 
gradually with height, roughly like the zonal wind con¬ 
tours (Fig. la). Moreover, the zonal flow on the giant 
planets is thought to have an approximately barotropic 
structure in the upper troposphere, so that critical lati¬ 
tudes of waves do not vary strongly with depth. Yet EMF 
appears to be peaked in the upper troposphere (Liu and 
Schneider 2014). Such arguments based on linear wave 
propagation, therefore, also do not seem to account for 
the entire EMF structure. Even if they would, the ques¬ 
tion would remain how the mean-flow structures come to 
be organized in such a way that EMF becomes concen¬ 
trated in the upper troposphere. The latter question also 
arises in the context of linear stochastic models that de¬ 
scribe fluctuations around a prescribed mean flow (Farrell 
and loannou 1996a,b). Such models are successful in re¬ 
producing midlatitude eddy statistics, including an upper- 
tropospheric EMF enhancement (Whitaker and Sardesh- 
mukh 1998; Zhang and Held 1999; DelSole 2001). But 
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the mean flow in these models is prescribed rather than 
influenced by the stochastic eddies, and characteristics of 
the stochastic forcing are fit, for example, to observations 
or GCMs. So it likewise is not clear how the upper- 
tropospheric EMF concentration arises. 

Despite their shortcomings, these appear to be the only 
two hypothesis that have been formulated to explain the 
upper-tropospheric EMF concentration: surface friction, 
and the greater linearity of the upper troposphere. In the 
present paper, we test these two hypotheses explicitly with 
the help of an idealized dry GCM, which captures the 
upper-tropospheric EMF concentration and allows us to 
investigate it systematically. The GCM is described in 
section 2. Section 3 tests the friction hypothesis by dis¬ 
cussing a circulation in which the poles are heated and the 
tropics are cooled, which exhibits EMF concentrated in 
the lower troposphere, although surface friction acts there. 
Section 4 tests the linearity hypothesis by comparing fully 
nonlinear simulations with quasi-linear (QL) simulations, 
in which nonlinear eddy-eddy interactions are suppressed. 
We show that nonlinear eddy saturation is more signiflcant 
in the upper than in the lower troposphere and that eddy- 
eddy interactions are essential for capturing eddy absorp¬ 
tion. In section 5, we propose a mechanism that accounts 
for the upper-tropospheric EMF concentration, focusing 
on the typical depth of baroclinic eddies, on the tropopause 
as a wave guide, and on nonlinear wave saturation in the 
upper troposphere. 

2. Idealized GCM 

The idealized dry GCM is based on the pseudo spectral 
dynamical core of the Geophysical Fluid Dynamics Fabo- 
ratory’s Flexible Modeling System. It is described in de¬ 
tail in Schneider and Walker (2006). The only modiflca- 
tion we adopt is the upgrade from the Robert-Asselin Alter 
to the Robert-As selin-Williams Alter in the leapfrog time 
stepping (Williams 2011). The primitive equations on a 
sphere are integrated using the pseudo spectral method on 
unevenly spaced vertical a-levels (Bourke 1974; Simmons 
and Burridge 1981). 

The model is thermally driven by Newtonian relaxation 
of temperatures toward the radiative equilibrium of a semi¬ 
gray atmosphere. The surface temperature T/ in radiative 
equilibrium is given as a function of latitude (j) by 

7/((/)) = 7/+A;,cos2<., (3) 

where Ah is the pole-to-equator temperature contrast and 
Tf the radiative-equilibrium surface temperature at the 
pole^ For Earth-like simulations, we set A/^ = 90 K 
and Tf = 260 K, with = Tf + 2A,,/3 = 320 K the cor- 
responding global-mean surface temperature in radiative 
equilibrium. The radiative-equilibrium skin temperature 
at the top of the atmosphere is set to Tf = 200 K. The 
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Fig. 1. Zonal wind (solid lines for eastward and dashed lines for 
westward winds, in ms“^) and (a) EMF (colors), (b) EKE (colors). 
The thick green line is the tropopause (2-Kkm“^ lapse rate). Based on 
ERA-40 reanalysis averaged from 1980 to 2001 (Uppala et al. 2005). 


radiative-equilibrium temperature decreases monoton- 
ically with altitude and is given as function of the optical 
depth (io(0) and the pressure p by 


r{pA) = Tl 


1 +do{(l)) ( — 
\Po 


0C-\ \ 


(4) 


where a = 3.5 is approximately the ratio of absorber (e.g., 
water vapor) scale height to density scale height. The op¬ 
tical depth is specifled as 


doW 



(5) 


to achieve a continuous monotonic decrease of radiative- 
equilibrium temperatures from the surface to a constant 
temperature at the top of the atmosphere. The Newtonian 
relaxation time scale T is a function of latitude (j) and ver¬ 
tical coordinate a, as in Held and Suarez (1994), 


t:{p,<I>) = ' + (Ts ' - '^i ‘)max ^0, cos^{<!>). 

^ (6) 

The constants T/ and are the relaxation times in the inte¬ 
rior of the atmosphere and at the surface in low latitudes. 
Earth-like simulations are carried out with Gh = 0.85, 
Zi = 50 days and = 10 days. The Newtonian relaxation 
scheme is presented in more detail in Schneider (2004). 
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The radiative-equilibrium temperature profile Ts is stat¬ 
ically unstable in the lower atmosphere. A convective 
parametrization redistributes enthalpy vertically and mim¬ 
ics the stabilizing effect of latent heat release in moist con¬ 
vection. When an atmospheric column is statically less 
stable than a prescribed lapse rate yPj, with dry adiabatic 
lapse rate = g/cp and 7 < 1 , its temperature is relaxed 
toward yPj on a time scale of 4 hours. Details can be 
found in Schneider and Walker (2006). The implied (but 
not explicit) latent heat release increases as 7 decreases; a 
value 7=1 corresponds to vertical entropy homogeniza¬ 
tion through dry convection. Earth-like simulations are 
performed with 7 = 0.7, corresponding to a lapse rate of 
6.9Kkm-'. 

Dissipation consists of hyperviscosity acting on tern- 
perature, vorticity, and divergence, and of momentum and 
dry static energy diffusion in a 2.5-km deep boundary 
layer (Smagorinsky et al. 1965). 

The simulations are performed at horizontal spectral 
resolution T85 with 30 vertical cr-levels. All time averages 
are performed over 600 days after 1400 days of spin-up. 

3. Heating the poles and cooling the tropics 

Latitude dependent radiative forcing on Earth intro¬ 
duces a vertical asymmetry of the troposphere because, 
at leading order, zonal-mean vertical shear is proportional 
to zonal-mean meridional temperature gradients (thermal 
wind blance). Here we examine a circulation in which the 
pole-to-equator surface temperature gradientJs reversed. 
This is achieved by setting A/^ = —90 K and T/ = 380 K; 
that is, radiative-equilibrium temperatures near the poles 
are larger than near the equator. The Newtonian relaxation 
time scale to this radiative-equilibrium temperature is uni¬ 
form with Ts = Ti = 40 days. All other GCM parameters 
are unchanged from the Earth-like simulation introduced 
in section 2 . 

The EME structure (Eig. 2a) is an upside-down version 
of the Earth-like structure (Eig. 2b). Maximal EME oc¬ 
curs inside the atmospheric boundary layer, where drag is 
acting on the flow. Surface drag is known to affect prop¬ 
erties of macroturbulent eddies and the general circula¬ 
tion, including EME amplitude, jet strength, and jet loca¬ 
tion (James and Gray 1986; Robinson 1997; Chen et al. 
2007; Liu and Schneider 2014). It has also been sug¬ 
gested that surface drag might at least partially explain 
the upper-tropospheric concentration of EME (Held 2000; 
Vallis 2006). Our GCM simulation indicates that friction 
alone is not responsible for the upper-tropospheric concen¬ 
tration of EME, because the maximal convergence occurs 
in the frictional boundary layer when temperature gradi¬ 
ents are reversed. 

Baroclinic eddies in the simulation are generated, as 
on Earth, in the extratropical troposphere, as can be di¬ 
agnosed from the time evolution of a simulation started 


from a slightly perturbed axisymmetric state with zonal 
winds and temperature structure equal to the zonal mean 
of the statistically steady state shown above. The Charney- 
Stern (1962) necessary condition for baroclinic instabil¬ 
ity is satisfied: although surface potential temperatures 
increase poleward, so that the surface temperature gradi¬ 
ent is reversed relative to Earth’s, the interior-tropospheric 
PV gradient is also reversed (negative) in midlatitudes. 
The quasigeostrophic potential vorticity (QGPV) gradient 
along isobars approximates the PV gradient along isen- 
tropes (Charney and Stern 1962) and is given by 

dyq = P-dyyU + fdp^^^^ . (7) 

The reversal of the QGPV gradient in the interior tropo¬ 
sphere arises because the stretching term (third term on 
the right-hand side) dominates the planetary vorticity gra¬ 
dient jS = dyf, where y = Rij) is the meridional distance 
coordinate. In other words, the slope of the isentropes, 
I = —dyO/dpO, flattens sufficiently rapidly with altitude 
in the interior troposphere that the stretching term reverses 
the sign of the QGPV gradient. In polar regions and in the 
subtropics, the QGPV gradient has the same sign as /3, so 
that the Charney-Stern criterion for baroclinic instability 
is satisfied only in midlatitudes. 

Eor baroclinic instability in a QG flow, the divergence of 
the zonal momentum flux is equal to a temporally grow¬ 
ing exponential times a weighted vertical integral of the 
QGPV meridional gradient, with a strictly positive weight 
(Held 1975). As a consequence, linear theory points to 
EME divergence where the QGPV gradient is positive 
throughout the column, that is, in our case, outside the 
baroclinic zone in midlatitudes. Angular momentum con¬ 
servation then implies EME convergence within the baro¬ 
clinic zone. 

Indeed, in our simulation with reversed temperature 
gradient, like on Earth, EME is converging in midlatitudes 
and diverging at low-latitudes (Eig. 2a); that is, EME is 
directed poleward between the subtropical dissipation re¬ 
gions and the baroclinic generation region. Because of 
surface drag, EME convergence in midlatitudes results in 
westerly surface winds, and divergence results in easterly 
surface winds closer to the equator. Westerlies and east¬ 
erlies are associated with the analogs of Earth’s Eerrel 
and Hadley cells, as indicated by the Eulerian streamfunc- 
tion in Eig. 2e (to be comapred with the Earth-like situ¬ 
ation in Eig. 2f). Eerrel cells are here thermally direct 
but eddy driven and very shallow because EME is con¬ 
centrated close to the surface. Hadley cells are thermally 
indirect and also eddy driven. Close to the surface, two 
thermally direct and thermally driven cells can be seen 
equator ward of 15°. 

Reversed meridional temperature gradients imply that 
eddy entropy fluxes, which are generally directed down 
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the potential temperature gradient (e.g., Kushner and Held 
1998; Held 1999), are in the mean equatorward. Because 
the vertical component of the EP flux (1) is proportional 
to the eddy potential temperature flux, it is directed down¬ 
ward. This is illustrated in Figs. 2a, b. It implies down¬ 
ward, rather than the usual upward, propagation of wave 
activity (Edmon Jr. et al. 1980). Wave activity then ac¬ 
cumulates near the surface in the extratropics, propagates 
equatorward and dissipates, leading to EMF and EKE 
maxima there, despite the drag (Fig. 2a,c). Wave activ¬ 
ity appears to propagate horizontally where farther ver¬ 
tical propagation is inhibited, be it by a solid boundary 
in the case of reversed temperature gradients or by the 
tropopause, an interface at which the static stability in¬ 
creases, in the Earth-like case. 

This simulation shows that upper-level EMF enhance¬ 
ment for an Earth-like simulation and the underlying 
asymmetry between the lower and the upper part of the 
troposphere cannot be attributed to surface friction alone. 
We have verified that this conclusion is robust when the 
strength of surface friction is varied. The latter has a 
strong effect on the near-surface zonal and eddy kinetic 
energies, but not on the EMF strength. Similarly, we have 
verified that the conclusions continue to hold when sur¬ 
face friction is only applied to the zonal-mean flow but not 
to the eddies. Both in Earth-like and reversed insolation 
simulations, the EMF structure is qualitatively unchanged 
when surface friction acts only on the mean flow, consis¬ 
tent with previous studies of the separate effects of friction 
on mean flow and eddies (Chen et al. 2007). In both cases, 
we observe that EMF strength is approximately doubled 
when surface friction acts only on the mean flow. This 
likewise shows that frictional damping of eddies in Earth¬ 
like situations is not responsible for the upper-tropospheric 
EMF enhancement. 

Our reversed insolation experiment was motivated by 
theoretical considerations. Nevertheless, it is worthwhile 
noting that it can describe the atmosphere of planets with 
high obliquity, which were recently studied by Ferreira 
et al. (2014). 


4. Quasi-linear simulations 

As explained earlier, it has been suggested that an upper 
troposphere that is more linear than the lower troposphere 
favors meridional propagation of eddies and may explain 
upper-tropospheric EMF enhancement (Held 2007). To 
test this hypothesis, we compare fully nonlinear to QL 
simulations for Earth-like parameters. A QL model only 
captures the linear behavior of eddies and their nonlinear 
interaction with a mean flow. Nonlinear eddy-eddy inter¬ 
actions are suppressed. 


a. Quasi-linear model 

Eddies are defined as departures from the average over 
longitude A at constant latitude (j) (zonal mean): 


Throughout this paper, the overbar denotes a zonal av¬ 
erage at constant pressure p and the prime the depar¬ 
ture from this average. In the G CM, we use a surface 
pres sure-weighted zonal average Ps{')/ps along o sur¬ 
faces because surface pressure acts as a density in a co¬ 
ordinates (whereas the flow is non-divergent in pressure 
coordinates). 

The QL approximation means keeping eddy-mean flow 
interactions and removing eddy-eddy interactions. The re¬ 
moval of eddy-eddy interactions in the dry GCM consists 
of modifying the momentum and thermodynamic equa¬ 
tions as described by O’Gorman and Schneider (2007). 
For example, the time tendency owing to meridional ad- 
vection of a field a = a + a' by a velocity v = v + v' can be 
decomposed as 


da 

dt 



_da _da' 
^ dy ^ dy 



— V 


fda' 
dy ' 


(9) 


The term v'da' jdy in (9) represents the advection of the 
eddies by the ed dies them selves. It can be decomposed 
into a mean part v'da' jdy and a fluctuating part. The re¬ 
moval of the eddy-eddy interactions consists of keeping 
only the former in (9): 


da da 
dt ^ dy 








( 10 ) 


As a consequence, all triad interactions involving only 
eddy quantities are removed, such that interactions be¬ 
tween the eddies and the mean flow are the only nonlin¬ 
ear interactions retained. This approach has received some 
attention since early studies of rotating large-scale flows. 
For example, small-amplitude wave activity conservation 
theorems have been derived within its scope (Charney and 
Drazin 1961; Eliassen and Palm 1961; Dickinson 1969; 
Boyd 1976; Andrews and McIntyre 1976, 1978). Keeping 
only eddy-mean flow interactions in the tendency equa¬ 
tion (10) is equivalent to linearizing the equation for the 
fiuctuating part a', keeping the equation for a unchanged. 
Hence the “quasi-linear” denomination for the resulting 
set of equations. However, it does not imply that all non- 
linearities are small because nonlinear eddy-mean flow in¬ 
teractions can be strong. 

The QL approximation is well-posed in the sense that 
it preserves inviscid invariants consistent with the order 
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Reversed insolation Earth-Like 



Latitude Latitude 


Fig. 2. Comparison of a circulation in which the poles are heated and the tropics cooled (left column) with an Earth-like simulation (right 
column). Colors show (a, b) EMF, (c, d) EKE, and (e, f) the mass streamfunction. The solid and dashed lines in (a-d) indicate westward and 
eastward zonal wind in ms“^, and the dark grey arrows the vertical component of the EP flux. Note that color scales are different in the reversed 
insolation and in the Earth-like conflgurations. The dotted lines in (e, f) represent three isentropes (295 K, 320 K, and 350 K) as an indicator of 
the thermal structure. The thick green line marks the tropopause (a 2-Kkm“^ lapse rate contour). The statistically identical northern and southern 
hemispheres are averaged, in addition to averaging zonally and temporally (over 600 days after a 1400 days spin-up to reach a statistically steady 
state). 


of the truncation: mass, angular momentum, entropy, en¬ 
ergy, entropy variance and, when applicable, potential en- 
strophy. First-order invariants are conserved because first 
moment equations are unchanged, and second-order in¬ 
variants are conserved because neglected third-order terms 
only redistribute second-order inviscid invariants among 
scales. 


In terms of a statistical closure, the moment or cumulant 
equations implied by the QL equations are closed at sec¬ 
ond order; third-order cumulants do not enter the second- 
order equations. So the QL approximation corresponds 
to truncating the hierarchy of cumulant equations at sec¬ 
ond order—an approximation that has recently been called 
CE2 and has been used to study the dynamics of barotropic 
jets (Marston et al. 2008; Tobias et al. 2011; Srinivasan 
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and Young 2012; Marston et al. 2015). Statistical Struc¬ 
tural Stability Theory (Farrell and loannou 1996a,b) and 
some kinetic theories of statistical physics (Bouchet et al. 
2013) are essentially equivalent to CE2. 

b. Mean zonal and meridional circulations 

O’Gorman and Schneider (2007) described the mean 
zonal circulation in the QL model and compared it with 
that in the corresponding fully nonlinear simulation. We 
extend here some of their findings to the EMF structure 
and the mean meridional circulation. 

Figure 3 shows the mean zonal wind and Eulerian 
mean meridional streamfunction. The QL model pro¬ 
duces upper-atmospheric jet streams above surface wester¬ 
lies (”eddy-driven jets”), associated with meridional cells 
extending from the surface to the tropopause. The cir¬ 
culation in the QL approximation is compressed in the 
meridional direction compared with the full model. The 
mid-latitude jet has a limited meridional extent and a sec¬ 
ondary jet appears at higher latitudes. This is consistent 
with the Eulerian streamfunction, which exhibits narrower 
Hadley and Ferrel cells. Moreover, weaker cells appear 
at higher latitudes, in association with the secondary jets, 
such that four cells occur in each hemisphere (the contour¬ 
ing in Fig. 3 does not reveal the weak high-latitude cells). 
O’Gorman and Schneider (2007) attributed the meridional 
compression of the meridional circulation to the fact that 
the suppressed eddy-eddy interactions isotropize the ed¬ 
dies and thus generally increase meridional scales (Stone 
1972; Shepherd 1987). 

Thus, the QL model successfully reproduces some fea¬ 
tures of the fully nonlinear general circulation, such as 
eddy-driven jets and meridional cells, together with a re¬ 
alistic thermal structure (three isentropes are indicated in 
Fig. 4). However, significant qualitative and quantitative 
differences suggest that the redistribution of angular mo¬ 
mentum by large-scale eddies is not realistic. We now ex¬ 
amine the structure of the EMF. 

c. Eddy momentum flux 

Averages of EMF in the full model and in the QL model 
are compared in Fig. 4. The top panel (Fig. 4a) depicts the 
well-known picture of momentum transported by eddies 
in the upper troposphere from the subtropics into midlat¬ 
itudes. Hence, EMF is converging above midlatitude sur¬ 
face westerlies and diverging above low-latitude surface 
easterlies. Idealized dry dynamics with a uniform surface 
reproduce the essence of the zonal flow and EMF structure 
in Earth’s troposphere (Fig. la). 

The bottom panel (Fig. 4b) shows that in the QL model, 
EMF exhibits a fundamentally different pattern. Consis¬ 
tent with surface westerlies around 30° and 75°, EMF is 
converging in the atmospheric column above. However, 
instead of a poleward EMF over much of the baroclinic 
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Fig. 3. Mass flux streamfunction (colors) and zonal wind in ms“k 
(a) Full model, (b) QL model. The thick green line marks the 
tropopause. As in Fig. 2, the flelds are averaged zonally, temporally 
and over both hemispheres. 


zone like in the full simulation, we observe in the QL sim¬ 
ulation equatorward and poleward transport on the pole- 
ward and equatorward flanks of the jet respectively. There 
are no well-defined upper-tropospheric extrema. Poleward 
EMF is maximal in the mid-troposphere, in sharp contrast 
to the full simulation. 

EKE, shown in Fig. 5, is maximal in the core of the 
jet, just below the tropopause, both in the QL and in the 
fully nonlinear model. Despite weak upper-tropospheric 
EMF, especially on the equatorward flank of the jet, EKE 
is still maximal near the tropopause in the QL simulation. 
Because the ratio of EMF to EKE can be interpreted as a 
correlation coefficient between meridional and zonal eddy 
velocity components, this implies that, contrary to the full 
simulation, u' and v' decorrelate in the QL upper tropo¬ 
sphere, especially equatorward of the main mid-latitude 
jet. Despite being of large amplitude, baroclinic eddies do 
not transport much angular momentum from the subtrop¬ 
ics to midlatitudes. This suggests that the shortcomings 
of the QL approximation in reproducing the EMF struc¬ 
ture do not arise because eddies would not reach the upper 
troposphere, for example, but because vertical propagation 
might not be captured adequately. Instead, the QL approx¬ 
imation does not capture the dissipation of eddies. To ob¬ 
tain a more precise picture of wave activity propagation 
and dissipation, we compute cross sections of the QG EP 
vector F (Eq. 1) and of the flux of QGPV q, which is pro- 
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Fig. 4. EMF (colors) and zonal wind (contours), (a) Full model, (b) 
QL model. Note the different color scales. The EMF amplitude in the 
full model is about twice as large as in the QL model. The dotted lines 
represent three isentropes (295 K, 320K, and 350K). The thick green 
line marks the tropopause (a 2-Kkm“^ lapse rate contour). As in Fig. 2, 
the fields are averaged zonally, temporally and over both hemispheres. 


portional to V • F through the Taylor identity (Edmon Jr. 
et al. 1980), 


Both F and v'q' are shown in Fig. 6. The meridional ex¬ 
tent of lower-tropospheric negative QGPV flux is compa¬ 
rable in the full and the QL models. The EP fluxes are 
qualitatively similar, roughly below the 295-K isentrope, 
where the QGPV flux is dominated by the vertical gradi¬ 
ent of the meridional eddy flux of potential temperature. 
Baroclinic growth and vertical propagation of wave activ¬ 
ity seem fairly well captured by QL dynamics in this part 
of the atmosphere. 

Betweeen the 320-K isentrope and the tropopause, the 
EP flux is very weak in the QL model. QGPV fluxes 
are mostly positive, whereas in the full model, signiflcant 
meridional EP flux occurs, and QGPV fluxes are negative, 
indicating absorption of eddies on the equatorward flank 
of the jet. The lack of absorption in the QL model, and 
even weak emission as suggested by the positive QGPV 
flux, accounts for the absence of enhanced EMF in the 
subtropical upper troposphere. It is consistent with large 
values of EKE being associated with weak EMF. 

Interestingly, the QL model performs better in the lower 
troposphere than in the upper troposphere. This appears to 


Fig. 5. EKE (colors) and zonal wind (contours), (a) Full model, (b) 
QL model. Note the different color scales. The EKE amplitude in the 
QL model is about twice as large as in the full model. See Fig. 4 for 
other details. 

contradict the view of a more linear upper troposphere and 
of a more turbulent lower troposphere (Held 2007). How¬ 
ever, considering time averages gives little information 
about the dynamical processes involved. To shed some 
light on the dynamics of the QL approximation and on the 
role of large-scale eddy-eddy interactions in eddy absorp¬ 
tion, we perform baroclinic wave lifecycle experiments. 

d. Lifecycle experiments 

A lifecycle experiment solves an initial value problem 
and aims at understanding the development and saturation 
of a growing disturbance in a baroclinically unstable zonal 
flow (Simmons and Hoskins 1978, 1980; Thorncroft et al. 
1991; Merlis and Schneider 2009). The initial condition 
we use here is a small-amplitude disturbance of normal¬ 
mode form with respect to the mean circulation of the full 
model (Fig. 3). We choose the zonal wavenumber ki = 6 
of the disturbance, corresponding to the fastest growing 
mode for the full model. Radiative and boundary layer 
parametrizations are disabled. Lifecycle experiments are 
run for the full and the QL models. In both cases, we use 
a normal mode of the statistically steady circulation of the 
fully nonlinear model. 

The initial disturbance is obtained with a version of the 
GCM linearized around the mean circulation of the full 
model. We follow a similar procedure as in Merlis and 
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Fig. 6. Quasigeostrophic EP vector (grey arrows) and QGPV flux 
(colors), (a) Full model, (b) QL model. The solid and dashed contours 
indicate zonal wind (as in Fig. 4). See Fig. 4 for other details. 


Schneider (2009). As an initial condition for this linear 
analysis, we perturb at all vertical levels the odd merid¬ 
ional spectral coefficients (from 3 to 83) of the vorticity 
field, corresponding to the zonal wavenumber ki. Only 
the spectral modes of all fields corresponding to ki are ad¬ 
vanced in time. Surface pressure is uniformly rescaled 
such that the amplitude A/? = {{ps — is reset to 

IPa when A/? exeeds 10Pa (brackets (•) indicate a global 
average and (•) a zonal average). Temperature and veloc¬ 
ity fields are rescaled accordingly using geostrophic bal¬ 
ance. The rescaling procedure is repeated a few times until 
the disturbance growth rate and phase speed remain con¬ 
stant. 

Energy cycle Before comparing wave activity propa¬ 
gation in the QL and fully nonlinear lifecycles, we briefly 
discuss the time evolution of EKE and of the two energy 
conversions involving EKE: the baroclinic conversion Ce 
from eddy available potential energy (EAPE) to EKE, and 
the barotropic conversion Ck from zonal kinetic energy 
(ZKE) to EKE (Lorenz 1955). The baroclinic conversion 
Ce is significant during baroclinic growth; the barotropic 
conversion Ck is significant during barotropic decay, cor¬ 
responding to zonal fiow acceleration through decaying 
eddies (Simmons and Hoskins 1978). The dominant term 
in Ck involves EME convergence (Lorenz 1955). 

The evolution of EKE as a function of time is shown 
in Eig. 7; Ce and Ck are shown in Eig. 8. The fully non¬ 


Time (days) 

Fig. 7. Global mean EKE as a function of time in the lifecycle ex¬ 
periments for the full model (solid black line) and the QL model (dashed 
red line). 


X 10“"^ 
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Fig. 8. Baroclinic (black) and barotropic (red) energy conversions 
Ce and Ck as a function of time in the lifecycle experiments for (a) the 
full model and (b) the QL model. 


linear simulation exhibits several cycles of growth, satu¬ 
ration, and decay of baroclinic eddies. The two first cy¬ 
cles, from day 0 to day 28 and from day 28 to day 38, 
correspond to what is discussed in Simmons and Hoskins 
(1978). However, the QL model does not exhibit as clearly 
defined lifecycles (Eigs. 7 and 8b); the time evolutions of 
both Ce and Ck are different. Eirst, there is cyclical large- 
amplitude conversion from ZKE to EKE, especially after 
day 45. (Small-amplitude conversion also happens in the 
full model, for example at day 28; see Eig. 8a.) Second, 
when Ck is negative, there is also relatively large conver- 












































10 


Journal of the Atmospheric Sciences 


sion from EKE to EAPE, which does not occur in the full 
model. 

We now look at wave activity propagation and absorp¬ 
tion during barotropic decay and at the termination of life- 
cycles, to understand why QL lifecycles are so different 
and to make the connection with the different EME struc¬ 
tures. 

Wave activity diagnostics We compute cross-sections 
of the EP flux and of the QGPV flux, as in the lifecycle 
studies of Thorncroft et al. (1991). Eor small-amplitude 
conservative eddies, the QGPV flux is proportional to the 
wave activity tendency, as can be seen from equations (2) 
and (11). In addition, in the WKBJ approximation, the 
EP flux is transporting wave activity at the group velocity. 
The applicability of the small-amplitude and WKBJ ap¬ 
proximations in the troposphere is questionable (e.g.. Pot¬ 
ter et al. 2013), Nevertheless, the QGPV flux can be in¬ 
terpreted as the tendency of a more general wave activity, 
defined for nonlinear eddies of any amplitude (Nakamura 
and Zhu 2010). 

When the barotropic conversion of EKE to ZKE is max¬ 
imal, the QL approximation captures the EP flux and wave 
activity tendencies fairly well (Pig. 9b). The QGPV flux in 
the middle to upper troposphere is negative, indicating the 
surf zone where eddies are absorbed. However, at the end 
of each lifecyle, when both barotropic and baroclinic con¬ 
versions are weak, there is only weak re-emission of wave 
activity from the surf zone in the full model (Pig. 10b, 
left panel), while there is very strong re-emission in the 
QL model, for example, at day 45 (Pig. 10b, right panel). 
Wave activity radiation from the surf zone for QL dynam¬ 
ics results in barotropic re-growth of eddies (positive Ck 
at day 45, Pig. 9a, right panel). Associated with positive 
barotropic conversion is baroclinic decay (negative Ce at 
day 45, Pig. 10a, right panel), consistent with an equator- 
ward eddy flux of potential temperature and a downward 
EP flux (Pig. 10b, right panel), unlike in the full model. 

The reason why wave activity on average is not ab¬ 
sorbed in the QL upper troposphere (section 4b) can 
now be better understood. Low-level baroclinic growth 
and vertical propagation of wave activity toward the 
tropopause are well captured. But after wave activity is ab¬ 
sorbed in the upper-tropospheric surf zone, it is re-emitted. 
EP fluxes radiating from the tropical upper troposphere are 
particularly clear (Pig. 10b, left panel). Similar absorp¬ 
tion followed by re-emission also occurs on the poleward 
fiank of the jet (not shown), although absorption there is 
of larger amplitude than emission and explains why EME 
is concentrated on this fiank of the jet in a statistically 
stationary state (Pig. 4b). Inefficient absorption of wave 
activity but similar baroclinicity compared with the full 
model explains why EKE is larger in the QL model than 


in the full model (section 4b): wave activity slushes merid¬ 
ionally in the upper troposphere, amplifying EKE and pos¬ 
sibly leading to resonances within an upper-tropospheric 
wave guide. 

The analysis of lifecycles confirms that lower- 
tropospheric dynamics are fairly well captured by the QL 
approximation and that the upper troposphere is more non¬ 
linear than the lower troposphere (section 4b). Bursts 
of baroclinic growth in the lower troposphere, consist¬ 
ing of alternate growth and decay of EP flux, are cap¬ 
tured. This contrasts to some extent with previous stud¬ 
ies of baroclinic lifecycles, in which low-level nonlinear 
saturation was invoked to explain why baroclinic con¬ 
version saturates (Simmons and Hoskins 1978; Held and 
Hoskins 1985; Thorncroft et al. 1991), as an essential 
step in a ’’saturation-propagation-staturation” paradigm 
for baroclinic wave lifecycles. Our lifecycle study sug¬ 
gests that QL mechanisms, like depletion of MAPE 
through baroclinic instability, might play a significant role 
in lower-tropospheric baroclinic growth saturation. How¬ 
ever, strongly nonlinear mechanisms appear essential in 
the upper troposphere. It remains to discuss why upper- 
level eddy absorption is not fully captured in the QL 
model. 

Potential vorticity rearrangement in the surf zone 
Eor adiabatic inviscid motion, potential vorticity is con¬ 
served on isentropes. Analyzing potential vorticity fields 
on isentropes can therefore be used to diagnose wave- 
mean flow interactions (Hoskins et al. 1985). Here we use 
the potential vorticity for the primitive-equations system 
(Rossby-Ertel, PV hereafter), rather than the QGPV. We 
evaluate it on the 350-K isentrope, which lies in the upper 
troposphere at low latitudes and in the lower stratosphere 
in polar regions (Pig. 3). 

In a barotropic framework, Rossby waves can be ab¬ 
sorbed quasi-linearly or nonlinearly (Held and Phillips 
1987). QL absorption occurs through the Orr mechanism, 
when a mean flow shears eddies in the same direction in 
which they are tilted and thereby transfers their energy to 
the mean flow (Parrell 1987; Lindzen 1988). As the ed¬ 
dies are sheared, the distance between vorticity extrema 
decreases, and vorticity anomalies become well approx¬ 
imated by the ratio of velocity anomalies to typical dis¬ 
tances between neighboring extrema. Since this distance 
shrinks in the shearing process, velocity anomalies, and 
with them EKE, have to decreases because of vorticity 
conservation, leading to EKE transfer to the mean flow. 
The non-conservation of EKE is a direct consequence of 
the linearized dynamical operator acting on eddy fields be¬ 
ing non-normal in the energy norm (Parrell 1987; DelSole 
2004). 

Eddy absorption, and the subsequent effect on the mean 
flow, is essentially a potential-vorticity mixing problem 
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Fig. 9. EP fluxes, QGPV flux, and PV when barotropic conversion of EKE to ZKE is maximal. Left column: full model; right column: QL 
model, (a) A reproduction of Eig. 8 with vertical dashed blue lines added to indicate when the EP flux and the QGPV flux in (b) are computed, 
(b) EP flux (grey arrows) and QGPV flux (colors). The EP flux and the QGPV flux are averaged over one day, between days 22 and 23 for the 
full model and between days 28 and 29 for the QL model. As in previous flgures, solid contours are for eastward winds and dashed contours for 
westward winds, with 10ms“^ increments and the bold line indicating the zero-contour. The green line marks the tropopause. (c) Corresponding 
Rossby-Ertel PV maps on the 350K isentrope. The arrows are for the isentropic-density weighted winds. 


(see Dritschel and McIntyre 2008 for a review), as ex¬ 
pressed in the Taylor identity (11). For linear decay, PV 
mixing is performed by the shearing in the zonal direction 
of meridionally propagating eddies. However, PV mix¬ 
ing can be nonlinear, and in general is so in planetary at¬ 
mospheres. Analytical theories were developed to under¬ 


stand absolute vorticity (or PV) rearrangement in nonlin¬ 
ear Rossby-wave breaking, notably the Stewartson-Warn- 
Warn (SWW) solution for barotropic inviscid critical lay¬ 
ers of small-amplitude waves in a constant-shear mean 
flow (Stewartson 1977; Warn and Warn 1978; Kill worth 
and McIntyre 1985). This theory predicts the formation 
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Fig. 10. As in Fig. 9 but at the end of a lifecyle, when the barotropic conversion of EKE to ZKE is minimal. 


of Kelvin cat’s eye structures consisting of closed stream¬ 
lines between fixed points near the critical lines, with a 
zonal wave number corresponding to the breaking wave. 
These structures are advecting PV anticy clonic ally at lead¬ 
ing order as a passive scalar, forming small-scale PV fil¬ 
aments that are rolling up around the center of the cat’s 
eye. Ultimately, this leads to vorticity mixing, which is 
here mediated by structures that are not zonally symmet¬ 
ric, in contrast to the linear decay. 


PV maps during maximal barotropic EKE-ZKE conver¬ 
sion are shown in Eig. 9c. They clearly show that wave 
activity absorption is nonlinear in the full modelThe 
point AI marks a thinning PV filament, as described for 
LCl wave breaking in Thorncroft et al. (1991). The re¬ 
gion near A 2 resembles the Kelvin cat’s eye of a critical 


^The increasing amplitude of PV extrema between Figs. 9 and 10 is 
worth noting and indicates that the numerical scheme does not conserve 
PV. The non-conservation is particularly striking in the QL approxima¬ 
tion, probably because PV is not transferred as efficiently to small scales 
where it can be dissipated by hyperviscosity. 
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layer, as predicted by SWW theory. Both phenomena in¬ 
volve eddy-eddy interactions and enstrophy cascading to¬ 
ward small scales. The filament A i constitutes the eastern 
flank of a structure also reminiscent of a cat’s eye. Roll¬ 
up of filaments is not visible on the 350-K isentrope we 
are showing, but it is visible at lower levels. Correspond¬ 
ing structures arise in the QL model (near B\ and B 2 on 
Fig. 9c). Barotropic linear theory of critical layers for 
small-amplitude waves indeed predicts the development 
of cat’s eyes (Dickinson 1970), before linear theory breaks 
down. QL dynamics cannot capture the subsequent thin¬ 
ning and roll-up of PV filaments, which are essential for 
vorticity mixing and the absorption of wave activity in crit¬ 
ical layers. However, QL dynamics does capture the for¬ 
mation of cat’s eyes (Haynes and McIntyre 1987). The 
baroclinic structures we observe here are more complex, 
but the simplified barotropic small-amplitude framework 
elucidates why they develop. 

In contrast to the full model, in the QL model, dipoles 
of positive and negative PV anomalies near B\ and B 2 
in Fig. 10 persist because eddy-eddy interactions are not 
available to excite zonal wavenumbers larger than that of 
the breaking wave to allow fllamentation. Haynes and 
McIntyre (1987) describe a similar phenomenon for the 
QL SWW solution. The rotation of PV anomalies around 
the center of the cat’s eye leads to alternate phases of ab¬ 
sorption and overreflection of the same amplitude (instead 
of decreasing amplitude implied by fllamentation in the 
fully nonlinear case). Re-emission of wave activity in 
our experiments is the baroclinic large-amplitude equiv¬ 
alent of the overreflection phase. When the positive vor¬ 
ticity anomaly is on the eastern side of the cat’s eye, the 
tilt is southwest to northeast such that EMF is poleward. 
When the vorticity anomaly is advected to the west of the 
cat’s eye, the tilt becomes southeast to northwest because 
it joins with the positive anomaly northwest of the cat’s 
eye. Hence, EMF is equatorward. This is the essence of 
the overreflection phase. Because the phase of absorption, 
in which the correlation between u' and v' is positive, is 
followed by a phase of overreflection, in which the cor¬ 
relation is negative, the correlation is on average close to 
zero, as observed in the statistically stationary circulation 
(section 4c). 

It is important to stress that nonlinear structures medi¬ 
ating PV rearrangement (cat’s eyes) have a large merid¬ 
ional extent, spanning much of the baroclinic zone in both 
the full and the QL simulations. The small-amplitude cal¬ 
culations of Dickinson (1970) and Haynes and McIntyre 
(1987) are consistent with our finite-amplitude computa¬ 
tion: similar cat’s eye structures arise in both the nonlin¬ 
ear and the QL simulations. The fundamental difference 
between the two cases lies in the details of vorticity dy¬ 
namics inside the cat’s eyes. 

PV maps at the end of a lifecycle are shown in Fig. 10c. 
Weak meridional gradients and small-amplitude zonal 


structures in the fully nonlinear simulation are a conse¬ 
quence of eddies having been absorbed. In the QL model, 
large-amplitude waves and strong PV gradients indicate 
wave resonance, resulting from wave activity re-emission. 
The absence of eddy-eddy interactions prevents PV flla¬ 
mentation and transfer of enstrophy and of wave activity 
toward small scales, where they can be dissipated. How¬ 
ever, the shortcomings of the QL model are already man¬ 
ifest within the framework of conservative dynamics. For 
barotropic SWW critical layers, this is shown in Haynes 
and McIntyre (1987). The exact role of diffusion in wave 
absorption in a more complex system is not clear. 

EMF is concentrated on the poleward flank of the jet 
in the QL simulation (Fig. 4). More work is required to 
understand this fact. We can conjecture that QL absorp¬ 
tion through the Orr mechanism is more efficient on the 
poleward flank of the jet. Ineed, PV maps do not show the 
formation of cat’s eye on this flank of the jet but only the 
shearing of eddies by the mean flow. 


e. The role of barotropic triads 

The lifecycle calculations show that wave decay is pri¬ 
marily nonlinear for Earth-like parameters. To determine 
what components of atmospheric turbulence are crucial, 
we have also considered a simplified GCM in which only 
barotropic triads are retained, while baroclinic-baroclinic 
triads and baroclinic-barotropic triads are neglected. Tur¬ 
bulence, here taken to mean transfer of inviscid quadratic 
invariants among scales, is contained in the barotropic 
mode only. Defining a vertical average with brackets [•] 
and a zonal average with an overline, the time tendency of 
a scalar a = d-\-a' due to the meridional advection by the 
velocity v = v + v' is integrated as follows (cf. Eq. 10): 


da _da _da' ,da ^da' 
dt ^ dy ^ dy ^ dy ^ dy 




( 12 ) 


Restoring the barotropic triads gives an EMF structure that 
is much closer to the full model. The upper-tropospheric 
enhancement, with realistic amplitudes, is captured, as can 
be seen in the upper panel of Fig. 11. This is consis¬ 
tent with EP vectors extending toward the subtropical up¬ 
per troposphere and showing eddy absorption at the equa¬ 
torward flank of the jet (Fig. 11b). Nevertheless, eddies 
are still compressed in the meridional direction, and sec¬ 
ondary eddy-driven jets are developing, though at signifi¬ 
cantly higher latitudes than in the QL model. Eddy absorp¬ 
tion on the poleward flank of the jet is reduced compared 
with the QL model but is still more significant than in the 
full model (not shown). Lifecycle experiments confirm 
that nonlinear saturation does occur on the equatorward 
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flank of the jet. Nevertheless, absorption on the equator- 
ward flank is less efficient than in the full model. Signifi¬ 
cant wave activity is still re-emitted and then absorbed on 
the poleward flank of the jet. As a consquence, EMF di¬ 
vergence is stronger than in the full model on the poleward 
flank of the jet. This is likely linked to the compression 
in the meridional direction of the circulation. Nonethe¬ 
less, it is clear that allowing the nonlinear interaction of 
barotropic triads suffices to obtain a much more realistic 
baroclinic EMF structure. 

It is not clear how the baroclinic waves interact with 
critical layers and why barotropic interactions suffice 
to give the baroclinic EMF structure. The baroclinic- 
baroclinic interactions are likely not important because 
vertical wavenumber one dominates the vertical struc¬ 
ture of the atmosphere: baroclinic-baroclinic interactions 
only can affect wavenumbers larger than two. How¬ 
ever baroclinic-barotropic triads are a priori important. 
Neglecting baroclinic-baroclinic triads and barotropic- 
barotropic triads, while retaining baroclinic-barotropic tri¬ 
ads, also captures upper-level enhancement. Thus, al¬ 
though dissipation of wave activity is not simply re¬ 
lated to the barotropic (vertically avergaed) flow, retain¬ 
ing barotropic interactions (barotropic-barotropic triads or 
baroclinic-barotropic triads) suffices to capture some eddy 
absorption. However, it does not capture the nonlinear dy¬ 
namics fully. 

5. Why eddy momentum fluxes are concentrated in the 

upper troposphere 

a. Summary and discussion 

In section 4, we have shown that EMF is not enhanced 
in the upper troposphere in the QL model with Earth-like 
parameters because wave activity is not absorbed on the 
equatorward flank of the jet stream. Baroclinic eddies can 
be absorbed through two fundamentally different mecha¬ 
nisms: nonlinear saturation and the QL Orr mechanism. 
Lifecycle experiments show that for Earth-like parame¬ 
ters, and for rapidly growing baroclinic eddies, nonlinear 
saturation is more relevant. While the Orr mechanism can 
be captured by QL dynamics, nonlinear saturation can¬ 
not. Essentially, eddy-eddy interactions allow wave activ¬ 
ity absorption through PV rearrangement in the surf zone 
and wave activity dissipation through enstrophy transfer to 
small scales. Without eddy-eddy interactions, wave activ¬ 
ity is primarily re-emitted from the surf zone, leading to 
excessive EKE in the upper troposphere. 

Yet it is possible to construct circulations in which eddy 
amplitudes or the geometry of the mean flow favor QL 
wave decay over nonlinear saturation, for example, by de¬ 
creasing the pole-to-equator contrast of the radiative forc¬ 
ing, by decreasing surface friction, or by making the planet 
rotate faster (not shown). Under such circumstances, the 



Latitude 


Fig. 11. Simplified GCM with wave-mean flow and barotropic 
eddy-eddy interactions, as described in section 4e. (a) Zonal wind and 
EMF (as in Fig. 4). (b) EP vector and QGPV flux (as in Pig. 6). 


QL approximation performs better and captures upper- 
level EMF enhancement more accurately. In some cases, 
then, fundamental aspects of the tropospheric EMF ver¬ 
tical structure can be understood in terms of QL dynam¬ 
ics, without having to consider nonlinear mechanisms for 
wave absorption. Wave activity is generated linearly in 
the lower troposphere and propagates vertically. Then, the 
tropopause acts as a turning surface for Rossby waves, 
trapping them in the troposphere and guiding them to 
propagate meridionally. This explains why no substantial 
EMF extends above the tropopause. In our experiment in 
which the poles were heated (section 3), waves are prop¬ 
agating downward, and the solid surface plays a similar 
role, resulting in EMF concentration near the surface. 

b. Depth of baroclinic eddies 

There is a vertical level below which EMF is weak 
(Fig. 4). The EP flux is primarily vertical in the lower 
troposphere because baroclinic growth dominates the dy¬ 
namics. One question arises: above which altitude does 
the EP flux acquire a substantial meridional component? 
That is, at which altitude does EMF become signiflcant? 
We suggest this level to be controlled by the typical verti¬ 
cal extent of baroclinic eddies (Held 1978; Schneider and 
Walker 2006). 
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To verify this hypothesis, we compare simulations in 
which the midlatitude tropopause height is set by con¬ 
vection with simulations in which it is set by the typi¬ 
cal depth of baroclinic eddies (Schneider 2004; Schnei¬ 
der and Walker 2006). This is achieved by changing the 
lapse rate Jgjcp to which the convective parametriza- 
tion is relaxing temperature profiles in our idealized GCM 
(section 2). As the convective lapse rate is reduced (7 
gets smaller), the tropopause rises, and convection be¬ 
comes increasingly important for the extratropical ther¬ 
mal stratification (Schneider and Walker 2006; Schneider 
and O’Gorman 2008; O’Gorman 2011), setting the height 
of the tropopause for 7 < 0.6. To resolve the upper tro¬ 
posphere and lower stratosphere well as the tropopause 
height increases, we perform these simulations with 60 
vertical a-levels, instead of 30 levels in the previous sim¬ 
ulations. 

The EMF for simulations with 0.4<7<0.9is shown 
in Fig. 12a-c, together with the FP vectors and the QGPV 
flux in dashed contours to indicate the depth of baroclinic 
eddies. The convective lapse rate parameter 7 has an im¬ 
portant effect on the FMF structure: as 7 is increased and 
baroclinic eddies become more important in controlling 
the extratropical thermal stratification, FMF become more 
peaked in the upper troposphere. 

For small convective lapse rates (e.g., 7 = 0.4), the EMF 
is strong over a large vertical extent and is maximal well 
below the tropopause (Fig. 12a). The EMF divergence 
(Fig. 12b) exhibits a particularly complex structure, in 
comparison to larger 7 , because it has two distinct local 
maxima: in the mid-troposphere (a 0.5) and near the 
tropopause (a ^ 0 . 1 ). 

EMF convergence also exhibits two corresponding 
maxima, but the one in the mid-troposphere is weak. To 
elucidate how this structure arises, we show in Fig. 12g the 
correlation coefficient between meridional and zonal ve¬ 
locity anomalies where the overbar here 

stands for a time and zonal average. It appears that 
the absolute value of the correlation coefficient only has 
one clear mid-latitude maximum, near the level of the 
mid-tropospheric EMF convergence/divergence extrema. 
Farge absolute values of the correlation coefficient indi¬ 
cate large EMF relative to eddy amplitude (EKE), and thus 
significant eddy absorption. Hence, comparing Figs. 12d 
and 12a suggests that the mid-tropospheric EMF conver¬ 
gence/divergence extrema are caused by the absorption of 
baroclinic eddies originating from lower levels, while the 
mechanisms responsible for the tropopause extrema are 
distinct. In fact, the EMF convergence and divergence 
near the tropopause correspond to northward and south¬ 
ward fluxes of PV (the potential temperature flux contribu¬ 
tion to the QGPV flux there is weak, unlike in lower tropo¬ 
spheric layers). The northward PV flux seems upgradient, 
and such fluxes have recently been argued to arise from 


nonlinear wave breaking (Birner et al. 2013). The cir¬ 
culation with 7 = 0.4 challenges this explanation because 
wave breaking regions appear well below the tropopause. 
Other processes may be responsible for the northward 
near-tropopause QGPV flux, such as local barotropic in¬ 
stability involving shallow modes. 

As the convective lapse rate ( 7 ) increases, the altitude 
where the absolute value of the correlation coefficient be¬ 
tween u' and v' in midlatitudes is maximal also increases, 
occurring at a = 0.55, 0.4, and 0.3 for 7 = 0.4, 0.7, and 
0.9 (Fig. 12g-i). This closely follows the deepening of 
baroclinic eddies, as shown by the tropospheric QGPV 
flux (e.g.. Fig. 12g-i). The level up to which the EP flux 
penetrates and the level of maximum correlation roughly 
coincide. The mid-tropospheric maximum of EMF diver¬ 
gence at the vertical level of maximum correlation disap¬ 
pears as 7 is increased from 0.4 to 0.7, leaving only the 
near-tropopause extrema. Apparently, the gap between the 
tropopause and the penetration depth of baroclinic eddies 
shrinks as 7 increases, both being equal for sufficiently 
large convective lapse rates ( 7 > 0.7, see Fig. 12a-c). But 
only for 7 = 0.9 does the vertical structure of the cor¬ 
relation coefficient resembles that of the EMF. The dy¬ 
namics that account for the near-tropopause EMF diver¬ 
gence/convergence extrema at 7 = 0.4, and contribute to it 
at larger 7 , do not leave a noticable signature in the cor¬ 
relation coefficient; they appear to be distinct from the 
lower-tropospheric baroclinic-eddy dynamics. A simu¬ 
lation with 7 = 0.4, in which the vertical resolution is 
halved, accurately captures the mid-tropospheric EMF di¬ 
vergence/convergence extrema but does not exhibit the 
near-tropopause maxima. This also points to distinct, shal¬ 
low mechanisms being responsible for the near-tropopause 
EMF divergence/convergence extrema for 7 < 0.7. 

To diagnose where the mean flow allows Rossby-wave 
propagation, and thus how deeply baroclinic wave can 
penetrate, we use refractive indices as defined in Seager 
et al. (2003): 


nl = 




kf 


u —3i{G)RCOS (j)/ki cos 0 


+ R^F{N) (13) 


where, in z-coordinates (Harnik and Findzen 2001): 


F{N)=f 


ez/2h 


dz 


e d 
V2 Tz 




(14) 


We denote with N and h the Brunt-Vaisala frequency and 
the pressure scale height. The refractive index is obtained 
by assuming a normal mode solution 
to the linearized QGPV equation around a zonal mean 
flow of QGPV q (A and ki are the longitude and zonal 
wavenumber). The resulting wave equation yields (Harnik 
and Findzen 2001): 

V^T^ + = forcing and disssipation. (15) 
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Fig. 12. (a-c) EMF (colors), QGPV flux (dashed contours for negative values and solid contours for positive values, in 10 ^ms in the 
northern hemisphere only), and the EP flux (grey arrows for midlatitudes in the southern hemisphere only), (d-f) EMF divergence (colors), QGPV 
flux (contours, as in (a-c)) and EP flux (arrows, as in (a-c)). (g-i) Correlation u' v' v'^) 2 , QGPV flux (contours, as in (a-c)) and EP flux (arrows, 
as in (a-c)). Note the change of scales for different values of 7 . (j-l) Rossby-wave refractive indices rir (colors) and zonal mean wind (contours, in 
ms“^). In all flgures, the thick green line marks the tropopause (a 2-Kkm“^ lapse rate contour). Columns from left to right: 7 = 0.4, 7 = 0.7 and 
7 = 0.9. 


The vertical coordinate in the two-dimensional laplacian 
operator of Eq. (15) is rescaled by the Prandtl ratio f/N. 
Positive values of indicate the possibility for wave prop¬ 
agation, whereas negative index preclude it. We compute 
the refractive indices for the fastest growing normal mode 
on the zonal- and time-mean circulation. 

Rossby-wave refractive indices are consistent with 
baroclinic eddies not reaching the tropopause for 7 = 0.4: 
the upper-tropospheric refractive index is negative near the 
jets with turning surfaces = 0 in the mid-troposphere 
(Fig. 12j). For 7 = 0.4, baroclinic eddies are shallower 
than the tropopause; they do not propagate sufficiently 


high vertically to reach it. This is consistent with the EMF 
maxima matching the depth of baroclinic eddies, and it 
further indicates that EMF near the tropopause is due to a 
shallow local mechanism 

To summarize, the concentration of EMF in the upper 
troposphere emerges for seemingly two distinct reasons. 
First, baroclinic eddies appear to be absorbed at a height 
that is determined by their typical depth. When this depth 
corresponds to the tropopause height ( 7 > 0.7), EMF is en¬ 
hanced in the upper troposphere. Second, near-tropopause 
dynamics are responsible for some meridional momen¬ 
tum flux (and even stronger EKE relatively to lower levels 
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baroclinic activity); the dynamical processes responsible 
for this are unclear but appear to be local to the tropopause. 


6. Conclusions 

Using an idealized dry GCM, we have investigated how 
the vertical structure of EMF is controlled and how its con¬ 
centration in the upper troposphere arises. In a simulation 
in which the poles were heated relative to the equator, we 
obtained an upside-down version of Earth’s tropospheric 
circulation with EME and EKE enhanced close to the sur¬ 
face. This shows that surface friction is not responsible for 
Earth’s weak EME near the surface, as had been suggested 
(Held 2000; Vallis 2006). Nonlinear lifecycle experiments 
are consistent with this conclusion in that they exhibit en¬ 
hanced EME in the upper troposphere, although surface 
friction in them is disabled. 

The upper atmosphere favors linear Rossby-wave prop¬ 
agation more than the lower troposphere. It has been sug¬ 
gested that this explains the EMF asymmetry between the 
upper and lower troposphere (Held 2007). To test this 
hypothesis, we compared a fully nonlinear model to a 
QL model, in which interactions between eddies and the 
mean flow are retained while nonlinear eddy-eddy interac¬ 
tions are suppressed. We have shown that the QL model, 
despite some success in capturing important aspects or 
planetary large-scale dynamics (O’Gorman and Schneider 
2007), does not reproduce the vertical EMF structure and 
its concentration in the upper troposphere. The reason is 
that eddy absorption in Earth-like parameter regimes is 
strongly nonlinear: nonlinear eddy-eddy interactions in 
the surf zone in the upper troposphere are important for 
the absorption of wave activity. Wave activity is re-emitted 
from the surf zone when eddy-eddy interactions are sup¬ 
pressed, as we saw in the QL baroclinic lifecycle exper¬ 
iments. This results in excessive eddy kinetic energy in 
the upper troposphere. That is, although the upper tro¬ 
posphere appears relatively linear as far as wave propaga¬ 
tion characteristics are concerned (Randel and Held 1991), 
it is more nonlinear than the lower troposphere with re¬ 
gards to eddy absorption. The QL model captures lower- 
tropospheric dynamics (e.g., baroclinic growth and its sat¬ 
uration) more faithfully than upper-tropospheric dynam¬ 
ics (e.g., the nonlinear surf zone). The QL model does 
not capture the essence of the upper-tropospheric dynam¬ 
ics because the surf zone in which nonlinear eddy-eddy 
interaction are important is not small in the meridional di¬ 
rection (as it is in some idealized flows like SWW); rather, 
it spans the meridional width of the baroclinic zone. An 
important consequence is that there is no clear separation 
between the latitude of generation of baroclinic eddies and 
the surf zone. EMF in this case occurs when dissipation 
occurs in the wings of eddies. 


To understand the vertical EMF structure, it is fruitful 
to think in terms of generation, propagation, and absorp¬ 
tion of wave activity. Wave activity is generated in the 
lower troposphere, propagates upward, then turns merid¬ 
ionally and is absorbed in the upper troposphere (Sim¬ 
mons and Hoskins 1978, 1980; Thorncroft et al. 1991), 
preferentially at a level that scales with the typical depth 
of baroclinic eddies. Wave activity propagation above the 
tropopause is inhibited because of the strong stabilization 
of the stratification above it, leading to the tropopause act¬ 
ing as a turning surface or wave guide (Thorncroft et al. 
1991). As a consequence, on Earth, the EMF structure is 
peaked in the upper troposphere because the tropopause 
height and the typical depth of baroclinic eddies coin¬ 
cide (Held 1982; Schneider 2004; Schneider and Walker 
2006; O’Gorman 2011). We have shown that when the 
tropopause is set by convection and baroclinic eddies do 
not reach the tropopause, EMF has a more complex struc¬ 
ture with a maximum well below the tropopause that cor¬ 
responds to the typical depth of baroclinic eddies. Ad¬ 
ditionally, there can be near-tropopause maxima of EMF 
divergence/convergence generated by distinct dynamics, 
whose origin is not entirely clear but which also play some 
role in upper-tropospheric EMF enhancement when the 
tropopause height is set by baroclinic eddies. Similarly 
but with reversed sign, in the simulation in which the poles 
were heated more strongly than the equator, wave activity 
is generated in the mid-latitude troposphere, propagates 
downward, then turns meridionally and is absorbed near 
the surface. In this case, the surface plays the role of a 
turning surface, with EMF reaching its maximum in the 
lower troposphere. 

A more complete understanding of the EMF structure 
in baroclinic atmospheres would require an explanation 
of why wave-activity absorption occurs preferentially at 
a height scaling with the typical depth of baroclinic ed¬ 
dies. Also, near-tropopause dynamics should be investi¬ 
gated further because they seem to be responsible for sig¬ 
nificant EMF and so affect the large-scale circulation. 

The results obtained with the QL GCM have implica¬ 
tions for the development of statistical closures for large- 
scale atmospheric dynamics. What favors quasi-linear or 
nonlinear decay, even in a simple barotropic framework, 
and what determines their relative importance when both 
can occur, is not fully understood. Unraveling this would 
be an essential step to develop such closures. Successes of 
second-order closures based on QL dynamics are likely 
explained by quasi-linear absorption being favored be¬ 
cause of the structure of the flow (Farrell and loannou 
1996a,b; Bouchet et al. 2013) or because of model tun¬ 
ing (Whitaker and Sardeshmukh 1998; Zhang and Held 
1999; DelSole 2001). The central role of eddy-eddy inter¬ 
actions in planetary macroturbulence does not contradict 
weakly nonlinear descriptions and theories, for example of 
the thermal structure (Randel and Held 1991; Held 2000; 
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Schneider and Walker 2006, 2008), in the sense that non¬ 
linear eddy-eddy interactions only are essential for eddy 
absorption. Somewhat ironically, the atmosphere looks 
more linear in the presence of nonlinear eddy-eddy in¬ 
teractions than in the QL approximation: eddies are of 
weaker amplitude with respect to the mean flow in the full 
model than in the QL model because their absorption is 
inhibited in the QL model. 
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